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TERRAIN  MESOROUGHNESS  DESCRIPTION  AND 
ITS  APPLICATION  TO  MOBILITY  AND  COVER 


Richard  A.  Weiss 
Mobility  Systems  Division 
U.  S.  Army  Engineer  Waterways  Experiment  Station 
Vicksburg,  Mississippi  39180 

ABSTRACT.  -^The  mobility  of  vehicles  traveling  off  the  road  has  always 
been  of  interest  to  military  planners.  In  fact,  the  survivability  of  military 
vehicles  on  a  battlefield  depends  to  a  large  extent  on  the  degree  of  mobility 
and  cover  afforded  by  the  terrain.  In  order  to  develop  accurate  measures  of 
the  mobility  and  cover  characteristics  of  terrain,  a  study  was  made  of  the 
statistical  description  of  terrain  elevation  variations  on  a  scale  smaller 
than  shown  on  a  topographic  map.  This  study  shows  that  a  useful  set  of  terrain 
descriptors  are  the  standard  deviations  of  the  detrended  elevation  and  its 
derivatives.  The  probabilities  of  a  vehicle  encountering  specified  values  of 
slope  or  finding  cover  in  a  specified  interval  are  calculated  in  terms  of  the 
terrain  roughness  descriptors.  These  probabilities  can  be  used  to  quantify 
terrain  areas  of  a  battlefield  and  develop  map  overlays  showing  patches  which 
indicate  the  degree  of  slopes  expected  to  be  encountered  and  the  amount  of 
cover  that  should  be  available.  Numerical  values  of  the  probabilities  are 
calculated  for  several  terrain  areas.  _ _ _ 


1 ■  INTRODUCTION .  A  significant  part  of  the  determination  of  the  sur¬ 
vivability  of  military  vehicles  on  a  battlefield  is  an  estimation  of  their  mo¬ 
bility  on  the  battlefield  terrain  and  an  estimation  of  the  degree  of  cover 
afforded  by  the  terrain.  The  detection  of  a  target  and  the  subsequent  fire¬ 
power  accuracy,  including  hit  and  kill  probabilities,  depend  on  the  percent  of 
a  target  that  is  exposed  to  fire  and  on  the  speed  at  which  a  vehicle  can  move 
across  the  battlefield  terrain.  This  paper  presents  a  method  of  estimating  the 
mobility  characteristics  and  available  cover  of  battlefield  terrain  and  speci¬ 
fies  the  terrain  parameters  required  to  accomplish  this. 


This  is  the  second  of  a  two-part  study  of  the  methods  used  for  a  quanti¬ 
tative  description  of  terrain  roughness.  The  first  part  considered  the  descrip¬ 
tion  of  terrain  microroughness  with  applications  to  the  prediction  of  the 
dynamic  response  of  military  vehicles  operating  on  rough  terrain. An  analy¬ 
tical  description  of  microroughness  is  necessary  for  the  design  of  track  and 
wheel  suspension  systems  and  also  for  the  design  of  optical  observation  and 
sighting  devices,  gun  stabilization  systems,  and  many  other  complex  weapons 
systems  that  are  part  of  modern  military  vehicles.  Many  of  the  analytical  pro¬ 
cedures  introduced  for  the  description  of  microroughness  can  be  used  to  describe 
large-scale  variations  of  terrain  elevation. 

The  second  part  of  the  terrain  roughness  study  deals  with  large-scale  ele¬ 
vation  variations  on  a  scale  appropriate  for  the  description  of  the  mobility,  D 

cover,  and  concealment  properties  of  a  terrain  area. 3  For  the  purposes  of  this  C 

paper  the  large-scale  terrain  elevations  can  be  separated  into  two  classes: 
mesoroughness  which  describes  the  terrain  elevation  variations  on  a  scale  between 
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microroughness  and  the  elevations  described  by  a  topographic  map,  and  macro¬ 
roughness  which  is  essentially  the  terrain  elevations  given  by  contours  on  a 
topographic  map.  The  mesoroughness  is  not  obtainable  from  a  standard  topo¬ 
graphic  map,  and  the  terrain  elevation  variations  of  the  mesoroughness  may  af¬ 
ford  cover  for  vehicles  on  a  battlefield. 

Aside  from  vegetation,  the  cover,  concealment,  and  mobility  characteris¬ 
tics  of  a  mesoroughness  elevation  profile,  as  seen  from  a  point  on  a  contour  of 
a  topographic  map,  are  determined  by  the  magnitude  and  character  of  its  eleva¬ 
tion  variations.  Therefore,  an  analytical  description  of  terrain  mesoroughness 
is  required  for  the  prediction  of  cover  and  mobility,  and  the  U.  S.  Army  Engi¬ 
neer  Waterways  Experiment  Station  (WES)  was  requested  to  develop  improved  analy¬ 
tical  representations  of  terrain  mesoroughness  and  to  incorporate  this  develop¬ 
ment  into  cover  and  mobility  indices  for  military  mapping  purposes. 

Measurement  of  Terrain  Elevation  Profiles.  A  study  of  terrain  roughness 
begins  with  the  measurement  of  an  elevation  profile.  A  simple  and  very  accurate 
method  of- obtaining  a  profile  uses  the  theodolite  and  surveyor’s  staff,  but  this 
method  is  very  time-consuming.  The  measuring  wheel  method  is  quicker  but  has 
the  disadvantages  of  disturbing  the  terrain  and  not  being  useful  in  very  rough 
terrain.  Terrain  elevation  measurements  done  by  WES  utilize  the  theodolite  and 
surveyor's  staff. 

Photogrammetric  methods  of  measuring  terrain  elevation  profiles  are  simpler 
and  faster  because  extended  areas  can  be  measured  by  areal  surveys. They 
have  the  further  advantages  of  not  affecting  terrain  conditions,  and  measuring 
all  terrain  features  such  as  vegetation,  roads,  rivers,  ditches,  etc.,  as  well 
as  elevations.  Photo gramme try  has  several  disadvantages  including  the  facts 
that  vegetation  often  impairs  terrain  elevation  measurements,  evaluating  photo¬ 
graphs  is  time-consuming  and  expensive,  and  the  resolution  of  elevations  is 
limited  especially  for  high  altitude  photographs. ^ 

Terrain  profiles  determined  by  any  method  are  presented  as  measurements 
of  elevation  at  discrete  points  along  some  predetermined  line.  The  effect  of 
this  discrete  scanning  is  to  remove  from  the  actual  elevation  profile  frequen¬ 
cies  higher  than  (AL)~1  where  AL  *  scanning  length  ,  so  that  scanning  can  be 
represented  as  a  low  pass  filter.  The  spatial  frequency  (AL)”1  is  essen¬ 
tially  the  Nyquist  frequency  associated  with  the  measurement  of  the  elevation 
profile  and  enters  the  calculation  of  power  spectra  through  the  spectral  window 
functions  for  the  slope  and  curvature.  Photogrammetric  methods  tend  to  under¬ 
estimate  the  values  of  the  standard  deviations  of  elevation,  slope,  and  curva¬ 
ture  and  their  associated  power  spectra  especially  in  the  high  frequency  regions. 

Macroterrain  Roughness.  A  battlefield  area  can  be  divided  into  a  grid  of 
square  areas  at  whose  vertices  the  terrain  elevations  are  specified.  These  ele¬ 
vations  are  recorded  on  a  topographic  map  at  uniform  spacings  of  generally  about 
10  or  100  m.  These  terrain  elevations  form  a  three-dimensional  surface  called 
the  macroterrain  roughness.  The  macroroughness  describes  the  large-scale  trend 
of  elevation  variations,  and  is  suitable  for  the  descriptions  of  cover  and  ve¬ 
hicle  mobility  over  large  areas. 7  Terrain  elevation  variations  on  a  scale  less 
than  the  100-m  grid  length  are  not  obtainable  from  a  topographic  map. 
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Mesoterrain  Roughness.  Superimposed  on  the  macroterrain  (the  trend  given 
by  the  elevations  specified  at  the  100-m  intervals  of  a  topographic  map)  are 
the  elevation  variations  on  a  scale  less  than  the  100-m  interval — the  terrain 
mesoroughness.  Troops  and  equipment  may  possibly  be  concealed  and  covered  with 
the  terrain  elevation  variations  of  the  100-m  interval  connecting  two  elevation 
contours  of  a  topographic  map,  whereas  a  straight  line  or  other  extrapolation 
of  elevation  between  these  two  points  would  not  suggest  this. 

The  slopes  occurring  within  the  mesoroughness  may  differ  considerably  from 
that  obtained  from  a  topographic  map,  and  therefore  the  actual  vehicle  mobility 
in  a  battlefield  area  may  be  considerably  different  from  the  predicted  by  the 
slopes  obtained  from  a  topographic  map.  Knowledge  of  the  mesoroughness  terrain 
elevations  will  be  critical  for  estimating  transit  time  of  vehicles  across  a 
battlefield.  Obstacles  of  various  sizes  and  shapes  are  also  expected  to  occur 
within  the  mesoroughness  and  these  will  affect  the  degree  of  mobility  in  a 
given  area.  The  mesoroughness  description  is  necessary  for  the  prediction  of 
cover,  concealment,  and  mobility  over  small  distance  scales  and  small  time 
scales. 

The  mesoroughness  elevation  varies  throughout  a  battlefield  area,  and, 
accordingly,  the  amount  of  cover  for  a  target  and  the  mobility  of  a  vehicle  can 
vary  rapidly  over  small  distances  (~10  m) .  Therefore,  it  will  be  of  advantage 
for  military  purposes  to  be  able  to  estimate  the  mesoroughness  for  the  areas  of 
a  battlefield.  A  number  of  mesoroughness  descriptors  need  to  be  developed  in 
order  to  do  this.  The  battlefield  is  divided  into  areas  within  which  the  meso¬ 
roughness  descriptors  have  essentially  the  same  values.  The  mesoroughness  de¬ 
scriptors  can  be  estimated  for  a  given  area  by  sampling  a  number  of  mesoterrain 
elevation  profiles  in  the  area. 

A  reduction  of  the  size  of  the  grid  spacing  of  the  macroterrain  roughness 
descriptions  would  adequately  incorporate  the  mesoroughness,  but  this  would 
require  large  computer  storage  capacity.  A  more  reasonable  procedure  is  to 
describe  the  mesoroughness  by  a  set  of  stochastic  variables. 

The  determination  of  the  descriptors  of  the  random  mesoroughness  compo¬ 
nent  requires  the  removal  of  the  trend  of  the  terrain  elevation  data.  A  de¬ 
trending  procedure  is  described  in  Part  II  of  this  paper.  With  a  properly 
selected  filter  constant  the  detrended  elevation  profile  gives  the  mesorough¬ 
ness  elevations  between  two  points  on  a  topographic  map.  The  detrending  is 
accomplished  by  a  computer  program  RFNWUD.3  It  is  assumed  that  the  mesorough¬ 
ness  irregularities  can  be  described  by  a  zero  mean  stationary  Gaussian  (normal) 
random  process. 8»9  The  mesoroughness  displacement  and  its  derivatives  will  be 
represented  by  zero-mean  Gaussian  processes  which  are  described  by  the  standard 
deviations  of  displacement,  slope,  curvature,  etc.  The  standard  deviations  of 
the  mesoroughness  displacement  and  its  derivatives  to  any  order  are  calculated 
from  a  detrended  elevation  program  by  the  computer  program  RFNWUD. 

Cover  and  Mobility  Characteristics.  Cover  refers  to  protection  from 
direct  weapons  fire,  while  the  term  concealment  refers  to  features  that  would 
interrupt  the  line  of  sight.  For  instance,  vegetation  may  afford  concealment 
but  not  necessarily  cover.  Vegetation  is  not  considered  in  this  report,  and 
only  cover  (concealment)  due  to  terrain  elevation  variations  is  treated.  The 
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degree  of  available  cover  afforded  a  target  is  measured  by  the  fraction  of  tar¬ 
get  height  that  is  covered,  i.e.  the  depth  of  the  hole  in  which  the  vehicle 
sits  compared  to  the  vehicle's  height. 

The  mobility  of  vehicles  depends  on  a  large  number  of  terrain  features  in¬ 
cluding:  vegetaion,  soil  strength,  terrain  slopes,  etc.  However,  the  primary 
terrain  roughness  variable  that  enters  mobility  calculations  is  the  slope  ex¬ 
pected  to  be  encountered  along  some  path. 10  Each  vehicle  is  associated  with  a 
critical  slope  which  it  cannot  negotiate. 

Because  the  mesoterrain  roughness  is  represented  by  stochastic  variables, 
the  average  mesoroughness  displacement  and  the  fraction  of  target  height  that  is 
covered  can  only  be  described  by  a  probability  theory.  This  is  also  true  of  the 
slopes  expected  to  be  encountered  on  a  battlefield  area.  With  the  assumption 
of  zero  mean  Gaussian  distributions  to  describe  the  mesoroughness,  the  calcula¬ 
tion  of  average  values  and  probabilities  of  encountering  specified  terrain 
characteristics  is  relatively  simple  and  is  done  in  Part  III.  These  average 
values  and  probabilities  are  expressed  in  terms  of  the  mesoroughness  descrip¬ 
tors  (standard  deviations  of  terrain  displacement  and  its  derivatives)  and 
can  serve  as  indices  to  delineate  patches  on  a  map  that  have  similar  cover  and 
mobility  characteristics. 

The  basic  objective  of  this  paper  is  the  development  of  new  mesoroughness 
terrain  descriptors.  The  objectives  and  scope  of  this  paper  are  shown  in 
Figures  la  and  lb. 


2.  STATISTICAL  DESCRIPTION  OF  MESOTERRAIN  ROUGHNESS.  This  section  de¬ 
velops  the  parameters  necessary  for  the  description  of  terrain  mesoroughness. 

It  is  assumed  that  the  mesoroughness  can  be  described  by  a  zero-mean  stationary 
Gaussian  random  process  which  is  obtained  from  a  measured  elevation  profile  by 
a  suitable  detrending  process.  The  frequency  content  of  a  stationary  random 
process  is  described  by  the  power  spectrum  or  alternatively  by  the  autocorrela¬ 
tion  or  autocovariance  functions. H  It  is  shown  that  that  autocovariance  func¬ 
tions  are  completely  determined  by  the  standard  deviations  of  the  mesoterrain 
elevation  and  all  its  derivatives.  Therefore,  the  basic  mesoroughness  descrip¬ 
tors  are  the  standard  deviations  of  the  mesoterrain  elevation  and  its 
derivatives. 

The  standard  deviations  of  mesoroughness  displacement  and  its  derivatives 
are  calculated  from  detrended  terrain  elevation  data.  The  accuracy  of  the 
values  of  the  numerical  derivatives  obtained  from  the  elevation  data  decreases 
for  the  higher  derivatives.  But  many  terms  (higher  derivatives)  are  required 
for  a  useful  power  series  expansion  of  the  autocorrelation  function  for  large 
argument.  Therefore,  the  power  series  method  of  calculating  the  autocorrelation 
function  and  the  power  spectrum  is  more  formal  than  practical,  and  for  numerical 
calculations  it  is  more  useful  to  represent  the  power  spectrum  as  a  polynomial 
whose  coefficients  are  determined  by  using  only  the  standard  deviations  of 
displacement,  slope,  and  curvature. 

Although  the  power  series  representations  for  the  autocorrelations  func¬ 
tions  are  only  of  formal  value  for  determining  the  power  spectra,  they  can  be 
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used  to  calculate  the  probabilities  of  finding  specified  mesoterrain  displace¬ 
ments  and  slopes  in  a  small  spatial  interval  of  the  terrain.  These  probabili¬ 
ties  are  important  for  the  description  of  cover  and  mobility  characteristics. 

Detrending.  Except  for  relatively  flat  areas,  a  macroterrain  elevation 
profile  cannot  be  represented  as  a  stationary  random  process  because  it  has  a 
trend,  i.e.  a  variation  of  the  statistical  parameters  along  the  length  of  the 
profile.  In  general  the  nonstationary  character  of  the  measured  elevation  pro¬ 
file  has  to  be  removed  by  a  detrending  procedure  in  order  to  obtain  the  standard 
deviations  of  the  resulting  me so roughness.  In  the  cases  where  no  trend  exists 
the  measured  elevation  data  can  be  processed  directly. 1  The  nonstationary 
property  of  an  elevation  profile  can  be  removed  by  a  detrending  procedure  which 
removes  the  long  wavelength  (the  trend)  components  as  shown  in  Figure  2a. 

The  mesoroughness  will  be  treated  as  a  stationary  random  process.  The 
mesoroughness  is  extracted  from  the  measured  elevation  profile  by  using  the 
following  exponentially  weighted  moving  average. 12 


S(x)  =  h(x) 


CO 

~  2 \J  [h(x  " 


a)  +  h(x  +  a) 


where 

C  =  mesoroughness  displacement 
h  ■  terrain  elevation 

X  *  detrending  constant  (filter  constant) 
x  =  horizontal  distance 
a  =  integration  variable 

The  detrending  procedure  given  in  Equation  1  removes  all  wavelength  components 
of  the  macroterrain  elevations  that  are  greater  than  X  .  A  computer  program 
RFNWUD  was  developed  to  accomplish  the  detrending. 3 

Plausibility  arguments  can  be  given  to  select  a  value  for  X  .  If  the 
measurement  interval  of  the  macroroughness  is  L  the  filter  constant  can  be 
taken  to  be 

X  -v.  L  (la) 


in  order  to  remove  the  trend  of  the  macroterrain  as  shown  in  Figure  2b.  Another 
possible  choice  for  X  would  account  for  the  elevation  difference  between  two 
points  on  the  macroterrain  of  a  topographic  map  as  follows 


where  H  -  elevation  difference  .  Equations  la  and  lb  represent  intuitive  pos¬ 
sibilities;  in  fact  the  values  of  the  parameter  X  should  be  chosen  to  produce 


agreement  between  the  probability  calculations  of  cover  and  mobility  given  in 
Part  III  and  the  measured  cover  and  mobility  characteristics  of  an  area. 


Mesoterrain  Roughness  Description.  Terrain  elevation  is  specified  at 
finite  intervals  of  horizontal  distance  and  is  generally  a  continuous  function 
whose  derivatives  are  discontinuous  at  the  points  of  measurement.  For  the  pur¬ 
pose  of  calculating  standard  deviations  the  derivatives  at  each  point  can  be 
assumed  to  be  the  slope  of  the  straight  line  segment  to  its  right  or  left  as 
shown  in  Figure  2c.  To  simplify  notation  the  n'th  derivative  of  the  mesoter¬ 
rain  displacement  ^(x^)  at  the  point  x  will  be  written  as 


Cn(xi} 


dng 

dx? 


(2) 


where  n  *  0,  1,  2,...  As  usual  the  zeroth  derivative  is  just  the  displacement 
itself 


W  =  C(xi} 


(3) 


The  root  mean  square  (rms)  values  of  these  derivatives  are  calculated  as 
follows 


N  N 


1=1  i=l 


N  NT 


[w]2  W 

[w]2  <5) 

[w]2  (6) 

[w]2  <7> 

[w]2  (8) 
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where 


Zq  «  rms  value  of  displacement 
Z^  *  rms  value  of  slope 

Z2  =  rms  value  of  second  derivative  (curvature) 

Z^  “  rms  value  of  third  derivative  of  displacement 
Zn  =  rms  value  of  n'th  derivative  of  displacement 
N  *  number  of  points  where  elevations  are  measured 

The  number  N  is  generally  large  being  of  the  order  of  hundreds  or  thousands. 

The  mean  values  of  the  macroterrain  displacement  and  its  derivatives  are 
calculated  as  follows 


n  w 

Si  N  Z-w  .  n  N  jLj  S/V 


(9) 


.  .  dx. 
i=l  1 


i=l 


The  standard  deviations  of  macroroughness  displacement  and  derivatives  are 
given  by^ 


o2  =  Z2  -  (cM) 

n  n  \  n/ 


(10) 


The  standard  deviations  can  also  be  written  as 
N 


°n" 


(ID 


i=l 


where 


W  =  Si(xi>  -  Sn 


(12) 


are  the  values  of  mesoterrain  displacement  and  its  derivatives  measured  from 
their  mean  values,  i.e.  their  values  for  a  zero  mean  process  since  =  0  . 
The  computer  program  RFNWUD  calculates  the  standard  deviations  of  mesoterrain 
displacement,  slope,  curvature,  etc.  from  measured  elevation  data  using  the 
formulas  1  through  10. 

Autocorrelation  and  Autocovariance  Functions.  The  autocorrelation  and 
autocovariance  functions  of  the  mesoterrain  elevation  and  its  derivatives  are 
defined  byll 
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N 


N  5n(xi)  Sn(xi 
i=l 


+  a) 


Cn(a) 


R  (a) 
n 


(13a) 


(13b) 


where 

R  (a)  *  autocorrelation  function  of  n'th  derivative  of  elevation 
n 

C  (a)  *  autocovariance  function  of  n’th  derivative  of  elevation 
n 

n  *  0,  1,  2,  3,  . 

The  autocovariance  functions  can  also  be  written  as^ 


N 


C 

n 


X  w  vxi 

i=l 


+  a) 


(14) 


where  ip  ( x )  is  given  by  Equation  12.  Because  of  the  stationarity  assumption, 
the  parameters  En  ,  on  ,  and  £**  ate  independent  of  the  interval  a  .  From 
Equations  8,  11,  13a,  and  14  it  follows  that 


R  (0)  m  Zn 

n  n 


(15a) 


cn(°) 


(15b) 


Useful  quantities  that  describe  the  random  variables  tp 
lation  coefficients^-^  n 


are  the  corre- 


rn(a) 


Cn(a) 

cn(°) 


(16) 


where  n  *  0,  1,  2,  3,  .  The  correlation  coefficients  describe  the  statis¬ 

tical  properties  of  the  mesoterrain  displacement  and  all  its  derivatives.  The 
correlation  coefficients  will  be  used  in  Part  III  to  calculate  the  probability 
of  encountering  specified  displacements  and  slopes  in  a  given  spatial  interval. 

The  autocorrelation  and  autocovariance  ?re  even  functions  and  therefore 
have  the  following  even  power  Taylor  series  expansions^ 
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Cn2)(0)  2  Cn4>(0)  4 

Cn(ci)  =  Cn(0)  +  ~ 2! -  a  +  ~T! - a 


(17b) 


where 


R^(0)  «  j  1  th  derivative  of  R  (a)  evaluated  at  a  =  0 
n  n 

C^^(O)  =  j '  th  derivative  of  C  (a)  evaluated  at  a  =  0 
n  n 


The  correlation  coefficients  are  given  by  Equations  16  and  17b  to  be 

1  Ci2>(0)  2  1  Cn4)(0)  4 

rn(°0  =  1  +  2T  c‘(0)  a  +  4l  C  “(0)  a  + 


(18) 


It  can  be  shown  that  the  autocovariance  function  and  its  derivatives  are  re¬ 
lated  to  the  standard  deviations  of  the  displacement  and  its  derivatives  in 
the  following  way3 
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expressions 

r  f  _ 

hold  for  the  R '  ' )  (0) 
n 

coefficients  with  the  o's  replaced 

Z '  s  . 


The  Taylor  series  of  the  autocovariance  functions  (Equation  17b)  and  the 
correlation  coefficients  (Equation  18)  can  therefore  be  written  as 


(24) 


(25) 


In  this  way  the  correlation  coefficients  can  be  evaluated  in  terms  of  the  stan¬ 
dard  deviations  of  mesoterrain  displacement  and  its  derivatives.  The  general 
expressions  for  these  functions  are 
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(26) 


(27) 

(28) 

(29) 

(30) 

(31) 


The  coefficients  appearing  in  the  power  series  expansion  of  the  correla¬ 
tion  coefficients  have  a  simple  physical  interpretation  if  it  is  noticed  that 
the  characteristic  wavelengths  for  the  mesoterrain  displacement  and  its  deriva¬ 
tives  are  given  as  follows 


X 

n 


(32) 


where 

X  ■  wavelength  of  n'th  derivative  of  the  mesoroughness  dis¬ 
placement.  Specifically, 

Xq  =  2iroo/o^  =  characteristic  wavelength  of  displacement 

X^  «  2tto^/o2  =  characteristic  wavelength  of  slope 

Xj  *  =  characteristic  wavelength  curvature 

X^  *  =  characteristic  wavelength  of  the  third  derivative 

and  so  on.  Therefore  from  Equations  31  and  32  it  follows  that 
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for  j  =  1 ,  2,  3, 

In  terms  of  these  wavelengths  the  correlation  coefficients  are  written 


as 


/  „  ,  1  / 2tt a\2  i  /2na  2ira  V 

rn(a)  =  1  '  ITl  — )  +  4T  Ix-T—J 

\  n  /  \  n  n+1/ 

1  / 2Tta 

'  Vx» 


2na  2na 


Xn+1  An+2> 


If  the  following  set  of  dimensionless  numbers  are  introduced 

T  -  2ffa 
n  \ 


(34) 
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the  correlation  coefficients  can  be  written  as 
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Written  out  in  full  the  expressions  for  the  correlation  coefficients  of 
mesoterrain  displacement,  slope,  and  curvature  are  respectively 
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with  the  restriction  that  Tj  <<  1  .  Figure  3  shows  typical  -correlation  coef¬ 
ficients  calculated  using  Equations  38  through  AO  and  the  computer  program  MESO. 
These  curves  show  how  the  calculation  breaks  down  for  large  spatial  distance 
(large  Tj)  because  only  a  finite  number  of  standard  deviations  are  accurately 
obtained  from  a  measured  elevation  profile.  In  other  words  the  use  of  a  trun¬ 
cated  series  restricts  these  results  to  small  values  of  spatial  distance  ct  . 

Power  Spectra.  The  power  spectrum  measures  the  frequency  content  of  a 
random  process.  It  is  defined  as  the  Fourier  transform  of  the  autocorrelation 

I  I  1  O 

function  as  follows-*--*-*-1-’ 
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where 


k  =  2ir!2 

n  =  spatial  frequency 
n  =  0,  1,  2,  3,  - 


The  definition  of  the  power  spectrum  is  physically  valid  only  if  the  average 
values  of  the  random  processes  have  physical  significance,  because  according  to 
Equation  13b,  the  average  values  appear  in  the  definition  of  the  autocor¬ 

relation  functions  R  (a)  . 

n 


However,  for  the  mesoroughness  description  the  autocovariance  functions 
Cn(a)  are  of  more  physical  interest  because  elevation  profiles  are  measured 
from  an  arbitrary  baseline  and  this  leads  to  arbitrary  average  values  ^  for 
the  mesoroughness  profile.  The  autocovariance  functions,  however,  describe  the 
mesoterrain  roughness  relative  to  the  average  values,  i.e.  it  takes  the  average 
values  to  be  the  baseline  and  therefore  describes  zero  mean  random  processes. 
For  terrain  roughness  descriptions  the  power  spectra  of  physical  interest  are 
defined  by  the  following  Fourier  transform  pairs 
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where 


k  =  2irft 


Equation  41a  is  the  formal  definition  of  the  power  spectra  for  mesoter- 
rain  displacement,  slope,  curvature,  etc.  This  equatic,  can  be  used  to  deter¬ 
mine  Pn(k)  provided  Cn(a)  is  completely  specified  for  all  a  ,  and  this  is 
never  the  case  for  actual  terrain  elevation  data.  In  particular  the  integral 
(Equation  41a)  converges  only  if  Cn(a)-»0  for  large  a  ,  i.e.,  Cn(a)  must  be 
known  for  large  a  and  must  approach  zero  in  this  limit.  This  precludes  using 
the  series  expansion  in  Equation  26,  which  is  valid  only  for  small  a  ,  because 
in  this  case 


p„<k>  •  v-  E  (-»\i  /“21e"toda 


where 


F.(k) 


00 

-p 


cos  (ka)da 


The  integrals  F j (k)  do  not  converge,  so  that  the  power  spectrum  cannot  be  ob¬ 
tained  from  a  power  series  expansion  of  the  autocovariance  function. 

The  use  of  a  finite  set  of  data  points  precludes  the  complete  determina¬ 
tion  of  Cn(a)  using  Equation  14  and  therefore  Pn00  cannot  be  determined 
for  the  full  range  of  frequencies  0  <  ft  <  00  by  using  Equation  41a.  A  finite 
set  of  data  points  results  from  the  finite  length  L  of  the  elevation  profile 
and  from  the  fact  that  the  elevation  profile  is  measured  at  intervals  AL  *  L/N 
so  that  the  limits  on  the  domain  of  definition  of  Cn(a)  are  AL  <  a  <  L  . 

Then  the  spatial  frequencies  have  the  following  bounds  L-1  <  ft  <  (AL)~^  ,  so 
that  the  lack  of  information  about  Cn  (a)  for  large  a  leads  to  a  lack  of  in¬ 
formation  about  Pn(ft)  for  small  SI  ,  while  the  lack  of  information  about 
(^(a)  for  small  a  leads  to  an  uncertainty  in  the  values  of  Pn(ft)  for  high 
frequencies.  The  upper  frequency  limit  (2AL)-^  often  occurs  in  information 
theory  and  is  called  the  Nyquist  frequency. ^ 

In  order  to  determine  the  behavior  of  Pn(ft)  for  ft  <  1/L  and  for 
ft  >  1/AL  a  model  approach  for  the  mesoterrain  roughness  power  spectrum  is 
adopted.  The  procedure  uses  the  following  mathematical  model  for  the  power 
spectra  of  mesoterrain  displacement,  slope,  curvature,  etc.l*2 


(45) 
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The  coefficients  C  ,  D  ,  and  E  of  the  three- 
parameter  power  spectrum  model  are  evaluated  from  the  values  of  aG  ,  °i  » 
and  02  that  are  obtained  from  detrended  (and  in  some  cases  undertrended)  ele¬ 
vation  profile  data  measured  at  intervals  of  length  AL  .  ^ ^  This  model  pre¬ 
dicts  the  five  basic  types  of  power  spectra  that  are  shown  in  Figure  4. 

Spectral  types  3  and  4  exhibit  no  trend,  so  that  for  these  cases  it  is  possible 
to  use  undertrended  elevation  data  to  determine  the  power  spectrum.  -*->2  The 
five  spectral  types  can  be  used  to  classify  mesoterrain  areas. 


3.  MESOTERRAIN  COVER  AND  MOBILITY  PROBABILITIES.  The  main  objective  of 
this  paper  is  the  development  of  a  quantitative  method  for  estimating  the  ef¬ 
fects  of  a  specified  mesoterrain  roughness  on  the  survivability  of  a  vehicle  on 
a  battlefield.  Survivability  depends  in  part  on  the  mobility  and  cover  charac¬ 
teristics  of  the  mesoroughness.  The  mobility  characteristics  are  described  by 
the  degree  of  slopes  expected  to  be  encountered,  while  the  cover  afforded  a 
vehicle  is  described  by  the  expected  amplitudes  and  widths  (wavelengths)  of  the 
hills  and  holes  of  the  mesoroughness. 

This  part  of  the  paper  uses  the  terrain  descriptors  defined  in  Part  II  to 
calculate  the  probabilities  for  encountering  specified  values,  and  ranges  of 
values,  of  mesoterrain  elevation  and  slope.  For  the  description  of  available 
cover  the  specified  value  of  the  mesoterrain  displacement  is  the  depth  of  a 
hole,  generally  equal  to  the  vehicle  height  or  larger,  which  will  afford  cover 
to  a  vehicle.  For  mobility  considerations  a  critical  slope  may  be  specified 
which  would  limit  a  vehicle's  performance. 

As  described  in  Part  II  the  mesoroughness  elevation  variation  and  its 
derivatives  are  described  by  stochastic  variables  whose  parameters  are  the 
standard  deviations  aQ  ,  ,  02  • . .  of  the  displacement  and  its  deriva¬ 

tives.  For  the  calculation  of  cover  and  slope  probabilities,  a  Gaussian  dis¬ 
tribution  is  assumed  for  the  random  variables  of  terrain  displacement  and  its 
derivatives. 

Two  kinds  of  probability  index  are  calculated  in  this  paper:  (a)  the 
probability  of  finding  a  given  mesoroughness  elevation  and  slope  at  a  fixed 
point  in  the  battlefield  area,  and  (b)  the  probability  of  encountering  a  speci¬ 
fied  mesoroughness  elevation  and  slope  in  a  specified  small  interval  of  the 
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battlefield  terrain.  The  first  probability  is  useful  in  the  case  where  the 
vehicle  motion  (or  travel  time)  is  not  important;  while  the  second  type  of 
probability  is  used  for  moving  vehicles  looking  for  a  covered  firing  position. 
The  two  situations  are  quite  distinct  and  an  area  having  large  probability  of 
cover  at  a  fixed  point  in  the  area  may  exhibit  a  relatively  small  probability 
of  finding  cover  in  a  specified  small  distance  interval. 

For  a  fixed  point  in  the  battlefield  the  probability  of  cover  can  be 
evaluated  if  the  standard  deviation  of  mesoroughness  elevation  0o  is  known 
for  the  area  which  contains  the  point.  The  probability  of  a  moving  vehicle 
finding  cover  in  a  specified  interval  will  also  depend  on  the  spatial  frequency 
content  of  the  mesoterrain  roughness.  Thus  a  relatively  smooth  area  (of  long 
wavelengths)  will  afford  less  cover  in  a  specified  interval  than  an  undulating 
area  (of  shorter  wavelengths)  with  the  same  standard  deviation  of  me  Thoroughness 
elevation.  The  calculation  of  the  probability  of  finding  a  specified  elevation 
(cover)  in  a  given  interval  will  be  shown  to  depend  on  as  well  as  aQ  ,  so 

that  for  the  case  of  a  moving  vehicle  looking  for  a  covered  firing  position  the 
standard  deviation  of  slope  is  a  critical  terrain  roughness  descriptor. 

The  probability  of  encountering  a  given  mesoterrain  elevation  or  range 
of  elevations  (cover)  in  a  specified  interval  of  distance  corresponds  to  finding 
cover  in  a  corresponding  specified  interval  of  time  given  by  at  =  pt  ,  where 
a  =  specified  distance,  p  =  vehicle  speed,  and  t  »  specified  time  .  If  the 
time  interval  is  taken  to  be  the  time  interval  between  successive  rounds  fired 
by  an  enemy  gun  ,  then  the  vehicle  had  better  be  able  to  find  cover  in  the 
distance  Of  ■  ptf  else  its  chances  for  survival  will  be  small.  Therefore  it 
is  of  practical  value  to  calculate  the  probability  of  finding  cover  in  a  speci¬ 
fied  distance  in  terms  of  the  roughness  descriptors  and  to  use  these  probabili¬ 
ties  as  an  index  to  delineate  areas  of  a  battlefield  having  different  degrees 
of  cover. 

In  a  similar  way  it  will  be  shown  that  the  probability  of  encountering  a 
specified  slope  in  a  given  interval  will  depend  on  the  standard  deviation  of 
the  second  derivative  02  as  well  as  on  the  standard  deviation  of  the  slope 
.  Therefore  for  mobility  problems  over  short  distances  in  a  battlefield  the 
curvature  roughness  parameter  is  as  important  as  the  slope  roughness  parameter. 
This  probability  is  important  to  assess  the  possible  values  of  slopes  that  may 
be  encountered  in  a  dash  for  cover  in  the  interval  0^  *  pt^  . 

Probability  Density  Functions  for  a  Point  in  the  Battlefield.  The  calcu¬ 
lation  of  probabilities  of  finding  specified  values  of  the  mesoterrain  eleva¬ 
tion  and  its  derivatives  at  a  point  in  the  battlefield  requires  the  calculation 
of  probability  density  functions  and  probability  distribution  functions. 

These  functions  are  commonly  used  in  probability  theory.  Of  particular  inter¬ 
est  to  military  problems  will  be  the  probability  of  encountering  specified 
values  of  elevation  and  slope. 

This  paper  assumes  that  the  mesoterrain  elevation,  slope,  curvature,  and 
all  higher  derivatives  are  independent  stationary  random  processes  whose  dis¬ 
tribution  about  their  mean  values  are  given  by  Gaussian  probability  density 
functions  defined  by 
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If  the  stochastic  variables  are  measured  from  their  mean  values  by 

•  £  -  £m  ,  the  zero  mean  Gaussian  distribution  is  given  by 

n  n  n 


PG(V 


(50) 


where  as  defined  in  Part  II  the  variables  £n  and  <|i  correspond  to  the 
detrended  elevation  profile  assumed  to  describe  the  mesoroughness  as  follows 

£  =  n'th  derivative  of  mesoterrain  roughness  elevation  measured  from 

arbitrary  level 

ij;  *  n'th  derivative  of  the  mesoterrain  elevation  measued  from  its  mean 
value 


Only  the  standard  deviations  o  are  required  to  describe  the  Gaussian  dis¬ 
tributions.  The  probability  density  functions  determine  the  probability  of 
encountering  a  specified  value  of  . 

The  probability  density  functions  have  the  following  properties  for 
n  *  0,  1,  2,  3,  ... 
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and  for  the  zero  mean  probability  density  functions 
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Typical  shapes  of  the  p_(£  )  and  p„(i|i  )  functions  are  given  in  Figure  5a. 

o  n  vj  n 

A  futher  descriptor  of  the  zero  mean  Gaussian  distribution  is  the  average 
of  positive  values  only  of  the  stochastic  variables 
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Therefore  on//  2ff  is  a  measure  of  the  average  amplitude  of  the  variation  of 
the  stochastic  variables  about  its  mean  value  of  zero. 

If  joint  probabilities  are  required  it  should  be  remembered  that  the 
derivatives  £  are  not  completely  arbitrary  in  the  sense  that  oD  ,  o-^  , 

02  ,  03  ,  etc.nare  not  completely  independent. Therefore  the  joint  distribu¬ 
tion  is  not  simply  the  product  of  two  or  more  proability  density  functions,  but 
requires  the  introduction  of  correlation  coefficients  between  the  various 
derivatives. ^ This  has  not  been  studied  in  this  paper. 

Probability  Distribution  Function  for  a  Point  in  the  Battlefield.  The 
probability  distribution  function  defined  for  the  zero  mean  Gaussian  distribu¬ 
tion  is 
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where  8  *  a  specified  value  of  $  for  n  *  0,  1,  2,  3, 

determines  the  probability  that  <^n<  6^  ,  i.e.,^^*^^ 
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The  probability  that  4n  >  0^  is  then  given  by 

Prob  (4.  >  0)  =  1  -  Pp(0  )  (59) 

n  n  G  n  ' 

The  probability  distribution  functions  for  Gaussian  distributions  are  often 
written  as 


W  =  erf 


where  the  error  function  erf(x) 


is  defined  as 
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A  graph  of  the  function  erf(x)  is  shown  in  Figure  5b. 

Probability  of  Cover  for  a  Vehicle  at  a  Point  in  the  Battlefied.  The 
zero  mean  Gaussian  probability  distribution  representing  the  mesoterrain  eleva¬ 
tion  variations  can  be  used  to  estimate  the  probability  of  finding  cover  for  a 
vehicle  at  some  point  in  a  battlefield.  The  condition  for  complete  cover  is 
assumed  to  be  that  the  vehicle  sits  in  a  depression  whose  depth  is  equal  or 
greater  than  the  height  of  the  vehicle  <  -  H  ,  where  H  -  vehicle  height 
(see  Figure  5c).  The  probability  for  finding  cover  is  v 

Prob  (*Q  <  -  Hv)  -  PG(-Hv) 


=  i  +  erf(-Hy/oo)  (62) 

=  \  -  erf(Hy/oo) 

A  typical  graph  of  the  function  PG(-HV)  is  given  in  Figure  5d.  The  relevant 
mesoroughness  parameter  is  the  standard  deviation  of  the  mesoterrain  elevation. 

Slope  Probability  at  a  Point  in  the  Battlefield.  For  slope  probabilities 
the  situation  is  somewhat  different  because  both  positive  and  negative  values  of 
these  parameters  must  be  considered.  The  probability  distribution  function  for 
slope  is  expressed  in  terms  of  a  critical  slope  <(;{  such  that  mobility  is  pos- 
sible  for  •  The  probability  for  slope  being  outside  this  range 

is 
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and  this  can  be  used  as  a  mobility  index.  This  function  appears  in  Figure  5e. 
A  similar  analysis  can  be  done  for  the  curvature  and  higher  derivatives. 
Wheeled  vehicles  have  »  0.3  and  tracked  vehicles  have  ^  *  0.45  . 

The  situation  of  a  moving  vehicle  seeking  cover  in  a  specified  distance 
and  of  encountering  a  specified  range  of  slopes  in  this  interval  can  also  be 
quantified  by  the  theory  of  probability.  The  probabilities  of  encountering 
specified  values  of  mesoterrain  elevation  and  slope  in  a  given  interval  can  be 
expressed  in  terms  of  the  autocovariance  functions  defined  in  Part  II  and  in 
terms  of  the  standard  deviations  of  the  elevations  and  its  derivatives. 

Probability  Density  Functions  in  an  Interval  of  Travel.  The  probability 
density  functions  for  encountering  specified  values  of  the  zero  mean  Gaussian 

random  processes  ,  given  by  i p  -  6  ,  in  a  small  interval  a  are  given 
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8  ,  in  a  small  interval  a  are  given 
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where  n  ■  0,  1,  2,  3,  ...refer  to  the  mesoterrain  elevation  and  its  deriva¬ 
tives.  These  probability  density  functions  can  be  rewritten  in  terms  of  the 
correlation  coefficients  defined  in  Equation  16  as  follows1 
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The  probability  density  functions  can  be  rewritten  in  terms  of  the  stan¬ 
dard  deviations  using  the  power  series  expansions  in  Equations  25  and  27  as 
follows 
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where  S  is  defined  in  Equation  31. 
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The  probability  density  functions  can  also  be  written  in  terms  of  the 
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dimensionless  wavelength  parameters  Tn  given  in  Equation  35  as  follows 
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where 
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For  reference,  the  first  few  4>  are 
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For  small  values  of  a  the  values  of  are  nearly  unity. 

The  probability  density  functions  (Equation  69)  for  an  interval  in  the 
mesoterrain  can  be  written  in  terms  of  the  Gaussian  probability  density  func¬ 
tion  as  follows 
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Consider  now  the  special  cases  of  mesoterrain  elevation,  slope,  and  curvature 

(75) 
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A  result  similar  to  Equation  75,  but  involving  the  power  spectrum,  has  already 
appeared  in  the  literature. From  Equation  75,  it  is  clear  that  in  addition 
to  the  elevation  Gaussian  probability  density  function  pq(60»o0)  ,  a  spatial 
frequency  term  /oQ  appears  in  the  expression  for  the  mesoterrain  elevation 
probability  density  function.  Therefore,  to  first  order  the  standard  devia¬ 
tions  oQ  and  0j  ,  determine  the  elevation  probability  density  function.  The 
standard  deviations  of  the  second  and  higher  order  derivatives  enter  to  a  les¬ 
ser  degree  through  the  function  <J>0  .  Likewise  Equation  76  shows  that  the 
probability  density  function  for  encountering  a  specified  slope  in  a  given  in¬ 
terval  is  proportional  to  o2/o^  in  addition  to  the  slope  Gaussian  probability 
density  function. 


Probability  Distribution  Function  in  an  Interval  of  Travel.  The  probabil¬ 
ity  distribution  functions  associated  with  the  probability  density  functions 
given  in  Equation  74  are 
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where  P„(0  )  is  related  to  the  error  function  as  in  Equation  60. 

g  n 

Cover  Available  in  a  Specified  Interval  of  Travel.  The  probability  of 
finding  cover  for  a  vehicle  of  high  Hv  in  a  specified  interval  of  the  meso- 
roughness  is  calculated  in  terms  of  the  probability  distribution  function 
(Equation  78)  as  follows 
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The  probability  for  finding  cover  in  a  specified  mesoterrain  interval  depends 
on  as  well  as  oD  through  the  ratio  o^/oQ  1/ XQ  .  In  other  words  Equa¬ 

tion  79  shows  that  for  the  same  values  of  oQ  the  mesoterrain  containing 
longer  elevation  wavelengths  affords  less  cover  in  a  specified  distance.  The 
standard  deviations  of  the  second  and  higher  elevation  derivatives  also  in¬ 
fluence  the  degree  of  cover  through  the  functions  defined  in  Equation  71. 


Degree  of  Mobility  Possible  in  a  Specified  Interval  of  Travel.  The 
problems  of  determining  a  mobility  probability  index  for  mapping  purposes 
amounts  to  calculating  the  probability  that  in  a  specified  interval, 

where  =  critical  slope  beyond  which  a  vehicle  cannot  go.  This  is  calcu¬ 
lated  as  follows 
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Then  the  probability  of  finding  slopes  less  than  the  critical  slope  is  given  by 

2ao.t’ 


Prob(l*il  <'"l)  ■  ■  T[  *1  ,rt(*iAi) 


(81) 


The  probability  of  encountering  a  critical  slope  in  a  mesoterrain  interval 
depends  on  02  as  well  as  on  through  the  characteristic  slope  wavelength 

02/ai  'v  1/XL  .  The  standard  deviations  of  the  third  and  higher  derivatives 
also  affect  the  calculation,  but  to  a  smaller  degree,  through  the  function 
defined  in  Equation  72. 


Average  Distance  Between  Covered  Positions  and  Between  Points  of  Critical 
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covered  positions  and  between 
value  occurs  can  be  estimated 
rain  elevation  wavelength  A0 
length  A^  that  are  given  in 
following  results 


two  adjacent  points  where  the  critical  slope 
from  the  values  of  the  characteristic  mesoter- 
and  the  characteristic  mesoterrain  slope  wave- 
Equation  32.  The  simplest  assumptions  give  the 
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where 

a  =  average  distance  between  two  adjacent  covered  positions  for  a  ve¬ 
hicle  of  height 

a^  =  average  distance  between  two  slope  mobility  failures  for  a  vehicle 
whose  critical  slope  is  . 

Numerical  Analysis  of  Terrain  Roughness  Probabilities.  Macroterrain  ele¬ 
vation  data  were  available  for  only  two  terrain  sites  -  Freiensteinau  and 
Wetzlar  both  located  in  West  Germany.  A  1-mile  section  of  macroterrain  eleva¬ 
tion  data  was  selected  from  each  site.  Smaller  sections  of  100-m  length  within 
the  one  mile  sections,  the  mesoterrain  elevation  profiles,  were  analyzed  and 
the  results  compared  the  descriptors  of  the  full  1-mile  section. 

The  standard  deviation  of  the  detrended  elevation  can  be  used  as  a  measure 
of  the  relative  roughness  of  the  Freiensteinau  and  Wetzlar  sites.  These  stan¬ 
dard  deviations  appear  in  Figure  6a  in  terms  of  the  reciprocal  of  the  detrend¬ 
ing  parameter.  In  terms  of  this  descriptor  the  Freiensteinau  site  is  about 
four  times  more  rough  than  the  Wetzlar  site.  Nevertheless,  as  seen  in  Fig¬ 
ures  6b  through  6d,  the  characteristic  wavelengths  A^  and  A^  are  roughly 
the  same  for  the  two  sites. 

The  power  spectra  of  the  terrain  displacement  were  calculated  using  Equa¬ 
tion  45  and  the  techniques  developed  in  References  1  and  2.  The  results  ap¬ 
pear  in  Figures  7a  and  7b  from  which  it  is  clear  that  several  different  types 
of  power  spectra  can  occur  for  the  100-m  sections  within  a  1-mile  terrain  ele¬ 
vation  profile.  Therefore  the  power  spectrum  type  will  vary  along  a  terrain 
section,  and  the  mesoterrain  roughness  power  spectra  is  expected  to  be  dif¬ 
ferent  from  the  macroterrain  roughness  power  spectra. 

Figure  8a  shows  the  Gaussian  cumulative  probability  for  terrain  eleva¬ 
tions  as  calculated  from  Equations  60  and  62.  Figure  8b  shows  the  results  of 
using  Equation  62  to  calculate  the  Gaussian  cumulative  probability  for  finding 
cover  for  a  vehicle  of  height  Ry  =  8  ft  in  terms  of  the  detrending  parameter 
A  .  The  detrending  parameter  enters  the  probability  calculation  through  the 
function  0o(A)  given  in  Figure  6a.  Equation  79  is  used  to  calculate  the 
cumulative  probability  for  finding  cover  in  a  1-ft  unit  interval,  and  the 
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results  are  shown  in  Figure  8c.  A  peak  occurs  in  the  Freiensteinau  probabil¬ 
ity  curve  in  Figure  8c  because  the  probability  varies  inversely  with  XQ  as  in 
Equation  79,  and  \Q  is  a  rapidly  decreasing  function  of  X--*-  as  given  in 
Figure  6b. 

Figures  8d  through  8f  give  the  probability  distribution  functions  for 
encountering  critical  slopes  using  Equations  63,  80,  and  81.  The  critical 
slopes  were  taken  to  be  0.3  for  wheeled  vehicles  and  0.95  for  track-laying 
vehicles.  The  standard  deviation  of  the  slopes  at  the  Wetzlar  and  Freinsteinau 
terrain  sites  are  about  an  order  of  magnitude  smaller  than  the  critical  slopes 
so  that  the  probability  of  encountering  these  critical  slope  values  is  vanish¬ 
ingly  small.  No  terrain  data  is  available  for  a  site  which  would  exhibit  high 
probabilities  for  encountering  the  critical  slopes.  In  order  to  obtain  some 
numerical  results  the  elevations  at  the  Freiensteinau  and  Wetzlar  sites  were 
arbitrarily  magnified  by  a  factor  of  10.  Figures  8d  through  8f  give  the 
critical  slope  probability  distribution  functions  for  this  artificial  situa¬ 
tion.  However  since  relative  probabilities  are  of  interest,  it  is  clear  that 
Freiensteinau  has  the  higher  probability  for  encountering  a  critical  slope 
value . 


Mobility  and  Cover  Map  Overlays.  The  cover  and  mobility  probabilities 
can  be  used  to  construct  military  map  overlays  on  which  areas  having  distinct 
mobility  or  cover  characteristics  are  isolated  to  form  a  patchwork.  Each 
patch  of  a  map  overlay  for  a  specified  vehicle  is  associated  with  a  probability 
index  of  cover  or  mobility  as  calculated  from  the  standard  deviations  of  the 
elevation  and  its  derivatives.  Several  elevation  profiles  are  measured  for 
each  area  to  determine  a  descriptive  set  of  values  for  oQ  ,  a\  ,  °2  >  and 

so  on.  In  this  way  map  overlays  can  be  produced  in  a  logical  fashion  from 
some  elevation  profiles  measured  in  each  area. 


4.  CONCLUSIONS.  This  paper  develops  a  formalism  for  developing  cover 
and  mobility  map  overlays.  The  formalism  is  based  on  a  rigorous  application 
of  probability  theory  to  the  description  of  random  terrain  elevation  data. 
Terrain  descriptors  obtained  from  elevation  profiles  measured  in  a  battlefield 
area  are  used  to  develop  cover  and  mobility  indices  for  the  purpose  of  de¬ 
termining  patches  on  a  map  overlay  having  distinct  cover  and  mobility 
characteristics . 

The  studies  of  mesoterrain  roughness  gave  the  following  conclusions: 

a.  The  mesoterrain  roughness  can  be  extracted  from  a  measured  elevation 
profile  by  a  detrending  procedure  with  a  proper  choice  of  filter 
constant  (Part  II) . 

The  basic  set  of  terrain  descriptors  required  for  the  complete  speci¬ 
fication  of  mesoroughness  are  the  standard  deviations  of  the  meso¬ 
terrain  elevation  and  its  derivatives  (Part  II). 

c.  The  probabilities  for  encountering  specified  ranges  of  mesoterrain 
elevation  or  slope  in  a  given  distance  can  be  determined  from  the 
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standard  deviations  of  mesoterrain  elevation  and  its  derivatives,  and 
these  probabilities  can  be  used  as  indices  for  determining  patches 
on  a  map  that  have  distinct  cover  and  mobility  characteristics 
(Part  III). 
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